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Abstract. We report ground-state calculations and Monte Carlo simulations 
for a lattice-gas model of the underpotential deposition of copper on Au(lll) in 
sulfatc-containing electrolytes. In a potential range of approximately 100~150 
mV, this system exhibits a (V3x V3) mixed phase with 2/3 monolayer (ML) 
copper and 1/3 ML sulfate. Our simulation results agree well with experimental 
results and with other theoretical work. 

1. Introduction 

Underpotential deposition (UPD) is a process whereby a monolayer (ML) or less 
of one metal is electrochemically adsorbed onto another in a range of electrode 
potentials more positive than those where bulk deposition would occur [1] . The 
UPD of copper on Au(lll) electrodes in sulfate-containing electrolytes has been 
intensively studied, both experimentally (see detailed discussion in Ref. [2]) and 
theoretically [3,4]. In cyclic voltammetry (CV) experiments the potential is 
scanned slowly while the current density is recorded. The most striking feature 
of the CV currents observed with a Au(lll) electrode in sulfate-containing 
electrolyte is the appearance of two peaks, separated by about 100~150 mV, 
upon the addition of Cu 2+ ions [5,6]. Typical CV profiles are shown in Fig. la 
[2], together with our simulation results. In the potential range between the 
peaks, the adsorbate layer is believed to have a (v^x \/3) structure consisting 
of 2/3 ML copper and 1/3 ML sulfate. This picture is based on a variety of 
both ex situ and in situ experimental evidence, as discussed in Ref. [2]. 

2. Model and theoretical methods 

Our lattice-gas model for UPD of copper on Au(lll) in sulfate-containing elec- 
trolyte is a refinement of the Huckaby-Blum model [3,4]. It is defined by a 
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standard three-state lattice-gas Hamiltonian [7,8], 
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Here nf e {0, 1} is the local occupation variable for species X [X=S (sulfate) 
or C (copper atom)], and the third adsorption state ("empty" or "solvated") 
corresponds to nf~nf=0. The sums Sa' an< ^ J2i run over au ^h- 

neighbor bonds, over all equilateral next-nearest-neighbor triangles [9], and 
over all adsorption sites, respectively, denotes the effective XY interaction 
through an Zth-neighbor bond, and runs over the interaction ranges. The 
change in electrochemical potential when one X particle is removed from the 
bulk solution and adsorbed on the surface is —fix- The sign convention is such 
that <I>xy>0 denotes an effective attraction, and /lx>0 denotes a tendency for 
adsorption in the absence of lateral interactions. The bonds that correspond 
to finite lateral interaction energies are shown in Fig. 2. For large separa- 
tions the interactions vanish, and <J>gg is an infinite repulsion corresponding 
to nearest-neighbor sulfate-sulfate exclusion ("hard hexagons" [3,4,10]). We 
emphasize that the <1>xy are effective interactions through several channels, 
including electron-, phonon-, and fluid- mediated mechanisms [7]. 

The electrochemical potentials in Eq. (1) are (in the weak-solution approx- 
imation and here given in molar units) related to the bulk concentrations [X] 
and the electrode potential E as 



Here zx (X = S,C) are the effective clectrovalences of sulfate and copper, R is 
the molar gas constant, T is the absolute temperature, and F is the Faraday 
constant. The quantities superscripted with a are reference values which con- 
tain the local binding energies to the surface. They are generally temperature 
dependent if the effects of rotational and vibrational modes are considered. 

The coverages of sulfate and copper are defined as 9x = N^ 1 ^ nf, where 
N is the total number of unit cells in the lattice. However, it is experimentally 
observed that sulfate remains adsorbed on top of the copper monolayer in 
the negative-potential region, rather than becoming reduced and entering the 
solution. A simple estimate of the sulfate coverage in this second layer can 
be obtained as 6g = adc(l/3 — 6$), which allows the difference between the 
first-layer coverage 6$ an d its saturation value of 1/3 to be transferred to the 



fix = A + RT\nj-±-zxFE. 
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top of the copper layer. The factor a is a phcnomenological constant expected 
to be between zero and one. Since the transfer of sulfate between the gold 
and copper surfaces does not involve an oxidation/reduction process, the total 
charge transport per unit cell during the adsorption/desorption process is q — 

(2) 

— e[z s (0s+9g ) + zc@c], where e is the elementary charge. In the absence of 
diffusion and double-layer effects and in the limit that the potential scan rate 
dE/dt^O [1], the voltametric current i per unit cell of the surface is the time 
derivative of q. Using differentiation by parts involving the relations between 
the E and ^x, Eq. (2), as well as the Maxwell relation d9 s /dp,c — d^c/dps, 
we find the current density i in terms of the lattice-gas response functions: 



which reduces to its standard form for a=0 [11]. 

Although the experimental studies are carried out at room temperature, 
the zero-temperature phase diagram serves as a quite accurate guide to the 
path in the (/Ss,Mc) plane that the isotherms should follow. At constant tem- 
perature and pH, two factors influence the path: the adsorbate concentrations 
in the electrolyte and the electrovalences. As seen from Eq. (2), p,$ and fie 
depend linearly on E, with slopes determined by z$ and zq, whose values 
must be determined from experiments. Here we use zc—+2, z$=— 2. Thirty- 
two ordered phases were found by applying the group-theoretical arguments of 
Landau and Lifshitz [12], nine of which (denoted by (Xx7)^ in Fig. 3) are 
realized as ground states for interactions in the region of experimental interest. 
The ground-state energies depend on fig and jic and the lateral interactions. 
The repulsive second-neighbor three-particle interaction <J> g g disfavors the pure 
sulfate (VJixv^f/ 3 phase, which has not been experimentally observed in this 
system. For a fixed set of interactions, the zero-temperature phase boundaries 
are exactly determined by pairwise equating the ground-state energies. In or- 
der to easily explore the effects of changing the interactions, a program was 
developed which numerically determines the zero-temperature phase diagram 
by scanning jig and jic and determining the phase of minimum energy [13]. 

To obtain adsorption isotherms and CV currents at room temperature, 
we performed Monte Carlo (MC) simulations on a 30x30 triangular lattice, 
using a heat-bath algorithm [14] with updates at randomly chosen sites. In 
order to avoid getting stuck in metastable configurations (a problem which is 
exacerbated by the nearest-neighbor sulfate-sulfate exclusion), we simultane- 
ously updated clusters consisting of two nearest-neighbor sites. Each data point 
was obtained as an average over 2xl0 5 Monte Carlo Steps per spin (MCSS), 
sampling at intervals of 50 MCSS and discarding the first 4000 MCSS. 
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3. Numerical results 



The zero-temperature phase diagram corresponding to the interactions used in 
this work is shown in Fig. 3. For large negative /2s only copper adsorption is 
possible, and the phase diagram is that of the lattice-gas model corresponding 
to the triangular-lattice antiferromagnet with next-nearest neighbor ferromag- 
netic interactions [15]. Similarly, in the limit of large positive /2s and large 
negative /2c the zero-temperature phase is the (\/3x\/3)(/ 3 sulfate phase char- 
acteristic of the hard-hexagon model [3,4,10]. The phase diagram for interme- 
diate electrochemical potentials is quite complicated. For its < —22 kJ/mol, 
no sulfate adsorption occurs in the first adlayer, while if /2c<— 18 kJ/mol, no 
copper is adsorbed. The (\/3 x -\/7)c/ 5 phase corresponds to experimental ob- 
servations in copper- free systems [16,17]. It is enhanced by the fourth-neighbor 
sulfate-sulfate attraction, $ss ' anc ^ * ne (V3xV3)l^ 3 phase is disfavored by the 
second-neighbor repulsive trios, <& S g. The (V3x V^)y^ mixed-phase region in 
the upper right-hand part of the diagram is relatively large, due to the nearest- 
neighbor attraction between copper and sulfate, 3> S c- The isothermal path is 
chosen such that the distance between points A and B (measured by the elec- 
trode potential) equals the peak separation in the CV current shown in Fig. la. 
The narrow strip of {V3x\/7)\^ phase lies closely above the (\/3x\/3)y^ phase 
and is quite sensitive to the fourth- neighbor attraction, $gg , and the second- 
neighbor trio sulfate repulsion, $gg. The nearest-neighbor copper repulsion 
causes the appearance of the (\/3x v/3)^g phase. 

The potential scan path corresponding to the CV current and adsorption 
isotherms shown in Fig. 1 is indicated by the dashed line in Fig. 3. With the 
aid of the ground-state diagram, it is easy to analyse the simulation results. 
Starting from the positive potential end (lower right in Fig. 3), where 6> c «0, 
8 S drops from its saturated hard-hexagon value of 1/3 to approximately 1/5 
in the (\/3 x \/7)c/ 5 phase region. This is in reasonable agreement with the 
experimental results for the same electrochemical system in the absence of 
copper, where the (\/3 x \/7)^ 5 structure has been observed, but not that of 

(V3x \/3)c/ 3 [16,17]. As E is scanned in the negative direction, the cations 
are attracted toward the electrode. Since there are strong effective attractions 
between the two types of adsorbed particles, the adsorption of copper induces 
the readsorption of sulfate: they both increase their coverages in a narrow 
potential range near 175 mV to form the mixed (\/3x\/3)2/3 phase. This is 
exactly the scenario proposed by Huckaby and Blum [3,4], and it is an example 
of the enhanced-adsorption phenomenon discussed by Rikvold and Deakin [7] . 
This phase remains stable in a potential range indicated by the separation of 
the two CV peaks, until the electrode potential is sufficiently negative that 



the sulfate adsorption on the gold surface is disfavored in comparison with 
completion of the copper monolayer. The replacement causes another sharp 
change in the surface coverages, corresponding to the left CV peak. However, 
there is experimental evidence that part of the sulfate desorbed from the gold 
surface is not reduced and dissolved, but rather remains adsorbed in a formally 

(2) 

neutral submonolayer on top of the monolayer of copper, with a coverage 8 S w 
0.2 [6]. This corresponds to a~0.6 in Eq. (3), which was used to obtain the 
simulated CV current and surface-charge densities shown in Fig. 1. Figure lb 
shows a comparison of the simulation results for the charge transfer and the 
integral of the experimental CV current. 

The agreement between the experimental and theoretical results is reason- 
able, except for large positive E, where the model predicts less copper and more 
sulfate on the surface than indicated by the experiments. The disagreement 
between the theoretical and experimental maximum currents may be due to 
defects on the electrodes used in the experiments. 
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Figure Captions 



Figure 1: Experimental (dashed curves) and simulated (solid curves) results, 
(a): CV currents for a positive-going potential scan at 2 mV/s. (b): Integrated 
charge density. Also shown are zq9q (long dashed), z s # s (dotted) and Zs(#S + 
4 2) ) (dotdash). 



Figure 2: The relative positions of copper (filled circles) and sulfate (triangles) 
corresponding to the effective interactions used in Eq. (1). The numbers are 
the corresponding values of <&x Y usec ^ m this work, given in kJ/mol. The 
interactions are invariant under symmetry operations on the lattice. 



Figure 3: The zero-temperature phase diagram. Solid lines arc phase bound- 
aries and the dotted line is the path along which the isotherms are calculated 
(positive E towards the lower right). The phases are indicated as (XxY)g s . 
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